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Abstract. We present a discrete inverse scattering transform for all ABS equations excluding Q4. The 
nonlinear partial difference equations presented in the ABS hierarchy represent a comprehensive class of 
scalar affine-linear lattice equations which possess the multidimensional consistency property. Due to this 
property it is natural to consider these equations living in an N-dimensional lattice, where the solutions 
depend on N distinct independent variables and associated parameters. The direct scattering procedure, 
which is one-dimensional, is carried out along a staircase within this multidimensional lattice. The solutions 
obtained are dependent on all N lattice variables and parameters. Wc further show that the soliton solutions 
derived from the Cauchy matrix approach are exactly the solutions obtained from rcflcctionlcss potentials, 
and we give a short discussion on inverse scattering solutions of some previously known lattice equations, 
such as the lattice KdV equation. 



1. Introduction 

The Inverse Scattering Transform (1ST) has been widely used as a mathematical tool for obtaining solu- 
tions of intcgrable nonlinear partial differential equations of two independent variables, such as the Korteweg- 
de Vrics equation (KdV), the sine-Gordon equation and the nonlinear Schrddingcr equation (see e.g. Ablowitz 
and Segur [3]). One uses the 1ST to solve the initial- value problem, where an integrable initial condition 
is usually given on the line t = 0, the Jost solutions are constructed, the time-dependent scattering data is 
computed, and the full time-dependent solution is then reconstructed. The 1ST naturally describes soliton 
solutions, where for the case of the KdV these correspond to rcflcctionlcss potentials. A rigorous discussion 
of the direct scattering of the time- independent Schrodinger equation along the full line can be seen in Dcift 
and Trubowitz [TU] and the case of the half- line in Fokas et al. |13j . 

The discrete 1ST dates back to the 1970's (to the best of the author's knowledge) with the works of Case 
and Kac [5] [5] , where the authors solved a direct discretisation of the Schrodinger equation on the half-line 
n > 0. Subsequently Flaschka [11] looked at the addition of continuous time evolution to the discrete 1ST 
in order to solve differential-difference equations arising from the Toda lattice. Other applications of the 
discrete 1ST to differential-difference equations can be seen in [37], [3] and [2J. A new scattering problem 
was derived in |29) from an "exact" discretisation of the Schrodinger equation. This problem has since been 
studied by various authors [7J |28j [17] for the cases of continuous and discrete time evolution, and as an 
eigenvalue problem for analytic difference operators. In 2010 a rigorous formulation of the discrete 1ST for 
this equation was given in [BJ, where the forward scattering problem was carried out on the line n = 
and new sufficient conditions were derived for the required properties of the transmission and reflection 
coefficients. 

In the present paper we build on the foundations of the machinery in [BJ to present a discrete 1ST for 
all ABS equations pQ, excluding Q4, in a complex multidimensional setting. The ABS list consists of nine 
partial difference equations with multidimensional consistency |26| . |19j . each of which representing a class 
of partial difference equations through point transformations of parameters and Mobius transformations of 
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the dependent variable. The list consists of H-type, Q-type and A-type equations. The H-type equations 

HI: {u -u){u - u) + p 2 - q 2 = (1.1a) 
H2 : (u-u)(u-u) + (p 2 -q 2 ){u + u + u + u)-p i + q A = (1.1b) 

/ 2 2 \ 

m 5 : Q(uu + uu) - P(uu + uu) + 5 ( P p( J j = (1.1c) 
where P 2 = a 2 - p 2 7 Q 2 = a 2 - q 2 
are all covered in this paper's scope. Here we have adopted the standard notation 

U . U nrni U . Un-\-l,mi " ^n,m+l; " ^n+l,m+l 

where n and m (s Z) are the independent variables of a two-dimensional lattice, with associated parameters 
p and g respectively. Of the four equations of Q-type we consider all except the topmost Q4, that is 

Qls : p(u — u)(u — u) — q(u — u)(u — u) + S 2 pq(p — q) = (1.2a) 
Q2 : p(u — u)(u — u) — q(u — u)(u — u) (1.2b) 
+ pq{p — q)(u + u + u + u) — pq(p — q)(p 2 — pq + q 2 ) = (l-2c) 

Q3.5 : P(uu + uu) - Q{uu + uu) - (p 2 - q 2 ) (^(uu + uu) + ^g^) = 0, (1.2d) 

where p = 'i a _ a 'i , Q = a1 a a and the parameters p := (p,P), q := (q,Q) lie on the Jacobi elliptic curve 
{(x, X) : X 2 = (x 2 — a 2 )(x 2 — b 2 )}. The two remaining A-type equations are related by straightforward 
gauge transformations to Ql$ and Q3 OJ and as such are not mentioned here as distinct examples. The choice 
of parametrisation that we have used was presented in |20j , in which each equation is parametrised by the 
common lattice parameters p and q, which occur naturally in all their solutions. This is advantageous for 
the degeneration procedure carried out in Section [3] 

Rather than simply considering the equations (jl.lj) and (|1.2j) as equations imposed on a two-dimensional 
lattice, the multidimensional consistency suggests that the most natural way to consider these equations is 
in an N-dimcnsional lattice with independent variables n\, and parameters p±„ p^. That is to 

say we require that each equation hold on every elementary quadrilateral throughout all N lattice directions, 
and thus the solution u of each equation will have dependence on each independent variable m, tt-n and 
each parameter p\, p^. This is the space in which we carry out the inverse scattering transform, and 
in Section [5] we rigorously develop the machinery for this, for the case of Q3,5. Here the initial condition is 
taken on an infinite one-dimensional staircase within this N-dimensional space. In Section [3] we perform the 
degeneration procedure between the Q-type and H-type equations, resulting in multidimensional solutions 
for all equations in the lists (|1.1[) and (|1.2[) . In Section [4] we consider the case of an arbitrary reflectionless 
potential, and in Section [5] we give explicit solutions to some previously known lattice equations, including 
the lattice KdV equation, which is not multidimensionally consistent. Finally we mention that in [4] it was 
explained how a well-posed initial- value problem for an arbitrary finite quad-graph can be immersed into a 
canonical initial- value problem for an N-dimcnsional unit cube. Thus the initial- value problem that we solve 
also solves that for an arbitrary two-dimensional quad-graph. 



2. Inverse Scattering Transform for Q3^ 

In this section we develop the scattering transform for Q3^, where we assume that the solution u and all 
parameters are complex. We consider the equation being imposed on each quadrilateral of an N-dimensional 
lattice, with independent variables n\, and associated parameters pi, Pat- For example if N = 3 

with independent variables n, m, I and parameters p, q, r respectively, then we impose that (jl.2dp hold in the 
(n,m)-plane with parameters p and q, in the (n, Z)-plane with parameters p and r and in the (m, Z)-plane 
with parameters q and r. This is permissible due to the multidimensional consistency of the equation, which 
is a well-known feature of partial difference equations, and has no analogue in the continuous case. 
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2.1. Linear Problem for Q3a. 

Consider a particular direction within the N-dimensional lattice, with independent variable and as- 
sociated parameter p k . Along this direction we then define the biquadratic [T] Jf associated with (|1.2d[) to 
be 

X Pk (u,u) := P k (u 2 +u 2 ) - (2p 2 k ~a 2 - b 2 )uu + A_, (2.1) 

where u denotes a shift of u in the rife-direction. This biquadratic plays a central role in the linearisation of 
(ll.2d[) , and in fact a convenient means of expressing (|2.1j) is by defining the function U through the first-order 
equation 

S 2 



UU = X Pk (u,u) = P k {u 2 +u 2 ) - (2p 2 k ~ a 2 ~ b 2 )uu + ^p-. (2.2) 



Here U is determined up to a multiplicative factor of a^ -1 ) k for some arbitrary constant a, and by imposing 
a = 1, U is determined uniquely. A Lax equation (first given in |21j ) in the n k direction for Q3,5 may then 
be written as 

P k u-{p 2 -b 2 )u C, 2 -b 2 

" 0. (2.3) 




'k 



UU 



Here C is the spectral parameter. Now since we could have chosen any one of the possible N directions we in 
fact have N biquadratics of the form (|2.1|) and N Lax equations of the form (|2.3p . The dependence of II on 
each of the N independent variables is defined through (|2.2|) . This in turn gives rise to tjN(N-I) Lax pairs, 

which are consistent (in the sense <p = <p) if and only if u satisfies (jl.2dp in each pair of lattice directions. 
2.2. Initial Staircase. 

We now consider the initial-value space for the direct scattering problem. Let T denote an infinite 
multidimensional staircase within the N-dimensional lattice, which is always non-decreasing in every lattice 
direction, and is defined through some stepping algorithm. An example for N=3 is a (1,2,4) staircase, that is 
one step in the n-direction followed by two steps in the m-direction followed by four steps in the ^-direction, 
and then repeating this process indefinitely. The staircase will then exist in a subspace (and span every 
direction of this subspace) of the N-dimensional lattice, and we denote this collection of lattice directions by 
I. 

Given the staircase L described above, we may enumerate the points along this staircase by a new inde- 
pendent variable i, and the parameter associated with each shift along L by pi — (pi. Pi). Thus each iteration 
i i + 1 will in fact correspond to an iteration in one of the N lattice variables, for example 713 — ^ + 1 
with parameter P3. We also use the notation <p to denote a shift of <f> along the staircase. 

By (|2.3[) we therefore have, along L, a sequence of Lax equations 

/ PiU-ipl-b^u C 2 -b 2 \ 

(p 2 -C 2 )^= [ _ U. (2.4) 

(Note that the parameters a and b, which appear originally in the elliptic parametrisation of Q3s, are constant 
for all lattice directions, and as such do not change along T.) This sequence of linear equations along T defines 
the direct scattering problem. Even though this initial-value space exists in up to N dimensions, the staircase 
itself, and thus the nature of the direct scattering problem, is one- dimensional. Allowing T to exist in multiple 
dimensions changes the parameter values along the staircase, but does not change the dimension of the direct 
scattering problem. As an analogy to the 1ST for continuous equations, the remaining N-l lattice variables 
play the role of N-l discrete "time" variables. 

Rather than consider the matrix problem (|2 .4[) it is convenient to consider the second-order equation 
satisfied by the first component of </>, which one can calculate to be 

'Pi+il - (p 2 +1 - p 2 )u - Piu\- 



(Pf + i ~ OHx ^= ~) 0i + ipi - OHi = 0. (2.5) 

The direct scattering problem for Q3^ is the solving of equation (|2.5[) . which will be discussed in due course. 
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2.3. Boundary Conditions for u and It. 

Soliton solutions to Q3s depending on an arbitrary number of independent variables and parameters are 
given in [50]. These are closely related to solutions of the NQC equation [53]. These solutions suggest that 
the most natural ansatz for the solution of Q3$ is 

u = AT (a, b)S A + T>T(a, -6)§ s + GT(-a, b)S e + 2XF(-a, -6)Sd. (2.6) 

where the exponential functions T are given by 

T( h\ TT ( (Pr + a )(Pr + b ) V"" / or . 



r=l 



and the four constants A, 23, C and D satisfy the single constraint 

S 2 

AD(a + b) 2 -'Be(a-b) 2 = — . (2.8) 

loab 

This final condition is necessary and sufficient for the background solution obtained by setting Syi — §s = 
Se = §d = 1 to solve Q3$. Thus from the ansatz (|2.6|) the problem of determining the unknown quantity 
u has been turned into that of determining the four unknown quantities S.a etc. We now impose some 
boundary conditions on these functions, namely that they satisfy 

$A ~ 1, §b ~ 1) §e ~ 1) §d ~ 1 as i — > — oo, 
S/i *~ S^, §2? ~ Sjj , §e ~ Sg, §d ~ Sjj as i — ^ +oo 

where S^t, S^, S e and S^, are constants satisfying S^S!^, = S^Sg = 1 (which ensures the background solution 
obtained by setting 8>a = S A , etc. does indeed solve Q3a.). This is the case for the soliton solutions [20j in 
which lower-order terms vanish exponentially, provided that we assume that the plane- wave factors appearing 
in the functions § A , etc. do not have modulus 1. 

With the above boundary conditions imposed on the functions § A etc., and therefore on u, we now 
calculate the boundary behaviour of the coefficient of <j> 1 in the direct scattering problem (12.51) . From the 
boundary behaviour of u, namely 

u ~ AF(o, b) + 3 T(a, -b) + &F(-o, b) + T>T{-a, -b) as i -> -oo, (2.9a) 
u ~ AF(o, 6)S^ + SJ( ffl , -6)S^ + eLF(-a, 6)S e + 2)J(-a, -6)S£, as i -> +oo, (2.9b) 

the boundary behaviour of It is then determined by (|2.2[) . As i — > — oo one may calculate this to be 

It - (a + b)AT(a, b) + (a - b)T>T(a, -b) - (a - &)£F(-a, b)-(a + b)DT{-a, -6), (2.10a) 

while as i — > +oo, 

It - (a + 6)AF(a, b)§>° A + (a - 6)SJ(a, -6)S£ - (a - b)£F(-a, 6)Sg - (a + 6)DJ"(-a, (2.10b) 

Therefore by combining (|2.9[) and (|2.10[) one finds that at both ends of the staircase the numerator of the 
coefficient of cj) 1 in (|2.5j) factorises and we have 

P i+ ll- (p? +1 -pf)u- PjU \ 

T = J ^ Pi + Pi+i as i — > ±00. (2-11) 

2.4. Forward Scattering Problem. 

The forward scattering problem for Q3a is the solving of equation (|2.5[) for the function <f>i. Due to 
the imposed boundary conditions on u (|2.9[) . which give rise to f|2 . 10[) and ()2.11j) . we make the following 
definition. 

Definition 2.1. The function v, which will be henceforth referred to as the potential, is defined to be 
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From (|2.11[) we clearly have v — > as i — > ±00. Equation (|2.5[) . which characterises the forward scattering 
problem for Q3§, then becomes 

(pUi - ( 2 )hi(i + 2; C) - (Pi+Pi+i + vt+i)Mi + 1; + fe 2 - C 2 )^i(*; = o. (2.13) 

The potential Vi is in fact the initial profile that must be specified before carrying out the forward scattering 
problem. As will be seen below this must satisfy a particular summability condition. Even though by f|2. 12[) 
this is equivalent to giving some nonlinear derivative of the actual solution u to Q3a, we shall see that, with 
the appropriate boundary conditions, it is sufficient initial data to be able to reconstruct the solution u 
through the inverse problem. 

2.5. Jost Solutions. 

We now look at particular solutions to equation (12.131) . 

Definition 2.2. The Jost solutions tp, ip, ip and which solve equation (|2.13p . are defined by the boundary 
conditions 



vM~nJi(E±§) 

where i n is some constant. 



> as % — y — oo, v y 



> as i -> +oo, (2-14) 



Since (|2 . 13[) is invariant under the map £ — > — C, it follows by the definition of the boundary conditions 
for the Jost solutions and uniqueness of the boundary value problem [T5], that 

0(i;O = ¥>(i;-C), ^(i;0 = ^(*;-0- ( 2 -15) 

Since the general solution to (|2.13[) involves two linearly independent solutions we may write 

tp = kip + Bip, 4' = kip + B(p (2-16) 

where A and B arc independent of i and A(£) = A(— £) and B(C) = B(— £). Moreover one can show that for 
any two solutions of (|2.13[) the Wronskian 

W{x,y) = { V l-C 2 ) 1 *{xy-xy) (2.17) 

is independent of i, and thus by choosing the relevant values at the boundary we have W(ip, ip) = W(ij), ip) — 
2C, implying 

AA — BB = 1. (2.18) 

2.6. Analyticity and Asymptotics. 

We now consider the analyticity and asymptotic properties of the Jost solutions and the spectral functions 
A and B. For convenience we make the following definition. 

Definition 2.3. The functions A, A and T, T are defined by 

-\n(^)^=\n(^) 4 

*- T n(!3) 4 - *-*n(l3) 4 - w 

For the purposes of this section it is also convenient to scale the potential v by introducing a multiplicative 
scaling factor A. That is, we rewrite the potential as 

v — > Xv, 

where A is a new complex parameter. This parameter is introduced solely for the analysis presented in this 
section, and thus will not actually appear in the solution u. Its benefit however is that it acts as a perturbation 
(from A = 1) of the initial- value space treated in the direct scattering problem. This perturbation does not 
affect the boundary conditions placed on u, and is introduced exclusively to analyse under what circumstances 
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a potential v will be admissible (see Definition I2.4[) . The Jost solutions and the spectral functions will (for 
the present) depend on A, and we make this explicit by writing A(i; Q — > A(i; A), etc. 

The linear equation (|2.13j) is a generalisation of that treated in [6J , where the direct scattering was carried 
out along a line rather than a staircase. This changes the parameters pi which are present along T but does 
not alter the nature of the scattering problem. Except for the adjustments required to allow for non-constant 
complex parameter values and the inclusion of the parameter A, the proofs of Lemma |2. II and Proposition 
2.21 follow those given in [BJ. 

Lemma 2.1. The functions A and T satisfy the following summation equation^: 

A(i; C, A) = 1 + A F & *•> OviHll C, A), T(i; £ A) = 1 + A £ G(i, I; C)^T(Z; £ A), (2.20) 

I— — oo l=i-\-i 

where 



F(i,l;t) 



1 

27 



1 



n 



Pr + C 



G(M;C) 



l 

2< 



n 



These summation equations have the Neumann series solutions 

A(i; £*) = ]£ A fe ^-(*; C). T(i; C, A) = ^ A fe J fc ( i; C) 



(2.21) 



fe=0 



k=0 



where 

H = l, 



H k+1 (i;C)= J2 F(i,l;QvtH h (!;Q, J = 1, J fe+1 (i;C) 



J=-c 



^ G(i, Z; ()viJ k (l; (). 

l=i+l 



Proposition 2.2. Let I denote the collection of lattice directions in which T exists, and assume that 

J2 N3ci i| (i + NI)<^ 



(2.22) 



w/iere 

% = 0C o (a, fo) > max 
Define the two open regions CP and CP* by 

- c 

r£l 



p r — a 



p r + a 



Pr + b 



Pr~b 



CP 



U{< 



Pr +C 



> 3C 



CP* := 



U{< 

rGl 



Pr + & 
Pr +C 



> 1. 



Pr - C 



then for C ^ ? and A e D = {A : |A- 1| < dj, d > 0, £/ie series solutions (|2.21[) converge absolutely. For 

each i the convergence is uniform in £ and A, and thus A and T are analytic in £ /or £ ^ {CPU9CP}, analytic 
in A /or X £ D and continuous in £ /or £ ^ CP. Equivalent results hold for A and T /or £ ^ CP*. 

The summability condition (|2.22[) is the only restriction that we impose on the initial profile, i.e. the 
potential v. This is a stronger condition that one finds for such partial differential equations as the KdV 
equation. The Q3^ equation represents a different class of equations to the KdV, and thus one should not 
expect that (|2.22p be equivalent to that obtained for the KdV equation [TU] or the lattice potential KdV 
equation [BJ . It is shown in Section [3] however that by the degeneration scheme given in |20j one can obtain 
solutions to all lower ABS equations, including the lattice potential KdV equation (HI) and the lattice KdV 



F(i, I; and G{i, I; Q) have removable singularities at ( = 0. To make this clear one may rewrite them in a different form, 
for example if T spans two lattice directions ni and U2 we have 



F(i,l;Q: 



Pi + C 



E 



Pi -c 
pi + c 



P2 +C 



E 



P2 - C 
P2+C 



/. 



for some t\ and satisfying t\ + ti, ■ 

2 Since %o > 1 both CP and CP* are bounded open regions with empty intersection. As an example, in the case of T being a line 
in the n-direction we have CP = {£ = x+iy : (x+\p re ) 2 + (y+\pi m ) 2 < (A 2 — l)|p| 2 }, CP* = {£ = x+iy : {x—\p re ) 2 + (y—\pim) 2 < 



(A 2 — l)|p| 2 }, where p ■ 



and A 



x'i-i' 
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equation. This degeneration process alters the summability condition (|2.22p . and for these equations this 
condition reduces to 

+00 

|ui|(l + |i|)<oo, 

i— — 00 

which is equivalent to that imposed for the KdV partial differential equation [TU] and the lattice potential 
KdV equation [BJ. For these reasons the condition (|2.22[) is not restrictive in the sense that it only allows 
the reproduction of the soliton solutions, but rather allows for a wider class of solutions with background 
radiation to be found. 



Corollary 2.3. 



X +00 . +OC 

A(C,A) = 1+- ]T vi?(l;C,\), B(C,A) = -- £ 

I — — OO 1 = — OO 



l-l 

n 



Pr - C 
Pr + C 



uiT(l;C,A). 



(2.23) 



A and A\ are continuous in C for C £ 7, except possibly at £ = 0, and analytic in C in the interior of this 
region. B and B\ are continuous in C for C £ {CP U CP*}, except possibly at £ = 0, and analytic in C m i/ie 
interior of this region. Both are analytic in A for A G D. 

Proof. (|2 . 23[) follow directly from comparing (|2 .20[) with 



i-l 



T-A+B I 



Pr+C 
Pr - C 



as 1 — > —00. 



Taking the Wronskian of (12.161) with respect to ip and (p gives 

A = ^W&, <P) = -£ [fa + C)T(i; C, A)A(i + 1; C, A) - fa - C)T(t + 1; C, A)A(i; C, A) 



B = — 



2C 



n 



Pr - C 

Pr + C 



A(i; C, A)T(i + 1; C, A) - A(» + 1; £ A)T(i; £ A) 



The remaining statements follow from these relations and Proposition 12.2 
Lemma 2.4. For A G D, 



A = 1 + - as |C| -> 00, T = 1 + O - as |C| -+ 00 



A = 1 + O I — ] as |d — > 00, B = O - as |C| -> 00 



(2.24a) 

(2.24b) 
□ 

(2.25) 
(2.26) 



Proof. From the series solutions (|2 . 2 1 1) it is straightforward to verify that Hk{i; C) = O ( <4j and ■/&(£; C) = 

O^^F^as|d— >oo, which proves (|2.25|) . To prove the results about A and B we use (|2.25[) in the expressions 
([233)1 . As |d -+ 00 we have 

l=1 ^ + £ v ,( 1+ o(r))=i + o(? 



I — — 00 



/ — — OO 



Z-l 

n 



Pr+C 
Pr - C 



«« 1 + =0 K , 



since for C £ {CPU CP*}, 



+ OO 

E 



i-i 

n 



Pr+C 
Pr - C 



'7 



< 2 |«,|3Cl'- 4 -l <oo. 

?= — OO 



□ 



Theorem 2.5. Assume that A does not vanish on 87 . Then either 

- A(C, 1) has a finite number of simple zeroes in C \ { CP U <9CP}, or 
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- There exists a punctured open disc D* around A = 1 such that for A £ D* , A(£, A) has a finite number 
of simple zeroes mC\{IPU 87}. 

Proof. A is analytic in this region, and since A ~ 1 + O f^J as |£| — > oo, one can show that no limiting 

sequence of zeroes exists and thus A must have a finite number of isolated zeroes in this region. For each 
such zero Cj there exists a locally convergent taylor series 

A(C,A) = 5>(A)(C-0(A))' 
i=i 

where each Q and each ai are analytic in A. If ai(l) 7^ then by continuity ai(A) 7^ in an open disc Dj 
around A = 1. If <xi(l) = then since ai(A) is analytic there again exists such a Dj on which di(A) 7^ for 
A 7^ 1. Letting D* be the intersection of all the sets Dj proves the result. □ 

Theorem 12.51 illustrates that the potentials which give rise to A(£, 1) having zeroes of order greater than 
unity in C\ {CPU90 3 } are exceptional cases, as the potential can simply be rescaled to avoid this circumstance. 
We therefore now ignore the parameter A altogether, and omit it from the notation, but make the following 
definition. 

Definition 2.4. Given a staircase Y, the potential v is called admissible if it satisfies 

+00 

]T 1^1^1(1 + \i\)< 00, 

i— — 00 

and gives rise to a spectral function A which has a finite number of simple zeroes mC\{TU 3CP}. 

Henceforth we assume that the potential has been scaled appropriately so that it is admissible. 

Corollary 2.6. Let Q be a zero of A in C \ {0 3 U 87}. If Q <£ 7* then ip(i; Q) = B(Q)(f(i; Cj)- If Q G IP* 
then ijj(i; Cj) — Cjip(i; Q) for some constant Cj. 

Proof. If Q i y* then the result follows from if) = Afi + Bip. If Cj € V* then by (|2.24ap we have W(ip, ip) = 0, 
which implies that these two functions are linearly dependent. □ 

2.7. Multidimensional Dependence of the Spectral Functions. 

We now consider how the Jost functions and the spectral functions A and B vary as functions of the lattice 
parameters n\, n^. By definition the functions A and B are independent of i, and by making the natural 
extension 

these functions are independent of every variable through which i cycles. For example if N=3 and we have 
a (1,1) staircase in the (n,m)-plane, then A and B will be independent of n and m. These functions may 
however have dependence on the remaining independent variable /. Consider one of the remaining N — I 
directions in which Y does not exist, with independent variable n s and parameter p s . From (|2.5p the evolution 
of 4>\ in this direction is governed by 

/ o 1 I / Pi(u — u) \ <? „ „ 1 

ivl - C 2 )Hi - I l ^ ; J <Pi + (pI - C 2 )Hi = 0, (2.27) 

where o denotes an iteration in the n s -direction. In accordance with the soliton solution behaviour |20j . we 
assume that in each of these remaining N-I lattice directions, 

P d - u) 

- 2p s as n s -> ±00, (2.28) 

U 

which implies that as n s — > ±00, 

(p 2 s - C 2 ) k k - tpsk + (pI - C 2 )Hi = 0. 
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The boundary conditions (|2.14[) defining the Jost functions however are not dependent on n s , and are thus 
inconsistent with this equation. 

Definition 2.5. The Jost functions which are consistent with (|2.27p in each of the remaining N — I lattice 
directions ( indexed by the set J ) are defined by 



,(N) 



Ps + C\ 2 S =(N) s TT (V* ~ C 



We now determine the dependence of A, B and the constants Cj on the remaining variables n s for all s G J. 
Proposition 2.7. The function A is independent of n s for all s G J: 

A(n.;C)=A(C). (2.29) 
The dependence of the function B and the constants Cj on n s , s G J, is given by 

B(n s ;C) = B(C) II > C ^ = <* U ^ 

Proof. For any given lattice variable n s , s G J, the functions <p( N ) and ip^ N ^ are linearly independent solutions 
of (|2.27p . so we may write 

Ps + C 



Ps-C 



</,( N ) = a ^ + /?<^ N) V = + II 

where a and (3 are independent of . Comparing this with ip = Atp + B^ shows that A = a is independent 
of n s , while the dependence of B on n s is a multiplicative factor of the form f ^ 5 ^ ^ . At ^ = Q with 
A((j) = 0, we have ip = Cjtp and thus Cj has the same dependence on n s as B. Carrying out this argument 
in all N — I directions gives the required result. □ 

Equations (|2.29[) and (|2.30|) are the main results of this subsection, and will be used in the following 
subsection where we carry out the inverse problem. 

2.8. Inverse Problem. 



We now consider the inverse problem, that is the reconstruction of the solution u to Q3s- In order to do this 
we first derive a singular integral equation for the N-dimensional Jost function A. Since we are considering 
this function as one depending not just on i, but on all lattice variables n\ to tin (i.e. A = A(ni, C))> we 
must utilise the multivariable dependence (|2.30p of B and the constants Cj. By (|2.29|) however, the function 
A is independent of all lattice variables. Furthermore since the scattering problem is symmetric in all lattice 
directions, there is no underlying "preference" given to any particular lattice variablq^- Moreover since the 
inverse problem is carried out in the complex plane of the spectral parameter £ only, it is convenient to alter 
the notation of the Jost functions by writing 

A(i;C)->A(0. 

Here it is understood that A still depends on all lattice variables m, hn, however for notational reasons 
we denote only the dependence on the spectral parameter 

Proposition 2.8. Let L be a closed contour in C \ {IP U 87} defining an interior bounded region 1Z + and 
an exterior unbounded region Hr . Assume that 7* C 1Z + , G TZ~ , that all the zeroes of A inC\{?U 87} 
are contained in 1Z + , and that all four points £ = ±a, ±6 lie in 1Z~ . Then given an admissible potential fj 



^this is not the case for partial differential equations such as the KdV equation, where there is a distinct difference between 
the dependence on x and t 
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on r, for £ 6 {7?. \ O 3 } i/ie N -dimensional Jost solution A(£) is i/ie unique solution to the singular integral 
equation 



a(0 = i-E 



i?(cr)A(a) 



(* + C) 



(C + (j 

where the sum is over the M zeroes of A in 1Z + , the constants Cj are given by 

C- ... _ „ B(C/ 



/o(ct) da. 



(2.31) 



if o £ r 



Ac(0) - ^ ~ " ' ~ J ' A C (0) 
and t/ie reflection coefficient and plane-wave factors are given respectively by 

N 



A(C)' 

Proof. Consider the relation (|2.16|) for A and T: 



Pr+C 
Pr - C 



T(Q 

A(C) 



A(C) - A(C) 



Pr + C 
C 



Pr 



n 



Pa +C 

c 



P« 



By choosing i D appropriately the plane-wave factor on the right may be rewritten as 



n 



Pr + C 
Pr ~ C 



n 



Ps + C 

c 



n 

rel 



Pr + C 
Pr~C 



n 



p« + C 
p*-C 



n 

r=l 



Pr + C 
Pr - C 



(2.32) 



and thus we have - A(Q = A(()R(()p((). The equation 



T(C) 
A(C) 



A(C) 
A* 



Si 



T(0 



k(Ci)(C-G) 



a/ 



A(0-i-E t 



i=i 



k(0)(C-0) 



A(C)i?(C)p(C), (2-33) 



then defines a jump condition between two sectionally holomorphic functions along the positively-oriented 
contour L. Given boundary conditions, the question of determining a function $(C) equal to these two 
functions on 1Z + and 1Z~ , which also satisfy the jump condition (|2.33[) on L, is a Riemann-Hilbcrt problem. 
Since the function defined on TV vanishes as |C| —> oo, the unique solution of this problem (see e.g. Gakhov 
[H] ch.II) is 

Using T(Qj) = CjA(Q)p(Q) or T(Q) = B((j)A(Cj)p(Cj) depending on whether or not Q e 7*, the solution 
$(C) for C £ TV gives (|2 .3 1[) . The existence and uniqueness of the solution to (|2.31j) follow from the existence 
and uniqueness properties of the corresponding discrete Gel'fand Lcvitan integral equation [6] that one may 
derive from this equation. In this scenario however we choose to use (|2.31|) rather than the Gel'fand-Lcvitan 
equation as the solution u of (jl.2d[) is most naturally expressible in terms of A. □ 

Remark 2.9. If all of the zeroes of A do not lie in !P* then consider a second positively-oriented contour L* 
lying inside TV such that the M zeroes of A lie between contours L* and L. By the residue theorem we have 

„<„),„ = £(, 



1 

2ttI 



R{a)A(a) 



which allows one to rewrite (|2.31[) as 

v(C) + v(0 



^ VA(G)(C + G) 



d\(a) = P (0, 



where the measure d\(a) is given by 2ird\(o~) = R(o~)da. By taking the continuum limit in which /o(Cfc) 



ikx-\-ik t 



, this integral equation becomes the linearisation of the KdV and Painleve II equations, first presented 



in [12) . This was used as the starting point for the direct linearization method |23j . and has since seen many 
generalisations and applications to other nonlinear partial differential and partial difference equations (see 
e.g. [H and ^). 
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We are now in a position to state the main result of this section, that is the reconstruction of the solution 
u of Q3^. The solution is expressed in terms of the Jost solution A which may be calculated from (|2.3ip . 

Theorem 2.10. Given an admissible initial condition Vi on the staircase T that is symmetric in its depen- 
dence on the parameters a and b, the N- dimensional solution u of Q3$ is given by 

u = AT(a, b)$ A + T>T(a, -6)S B + &F(-a, fe)S e + DJ(-a, -&)§©, (2.34) 

where J 7 (a, b) is defined by (|2 . T[) and 

§A=S(a,b), §s = §(o, — b), §e = S(— a, b), §d = §(— a, — b) (2.35) 

with 

§(a, b) = -^— b [{p k + a)A(o)A(6) - (p fc + 6)A(a)A(6)) , (2.36) 

where ^ denotes an iteration in any one of the N lattice directions, and and pk are the independent 
variable and associated parameter in this direction. 



Proof. Consider the two components of any one of the N Lax equations of the form ()2.4|) : 



4P fc u(c 2 -& 2 )y ri v u 

The equation governing <pi in the nfc-direction is 



(Pi - C 2 )^i - U) j 0i + (p£ - OHi = 0. (2.38) 

Now consider the points £ = ±6. At these points the equation governing (f> 2 in the n^-direction is 

(pi a 2 )^ 2 ( ^-""M & + (pg - a a )*fc = 0, 

thus at £ = b we may write 

(C 2 -b 2 )v { 2 N \C)} =c 1 ^ N \a)+c 2 ^ N \-a) 

(C - 6 ) Va (C) t h = cMa) 11 f— + c*p(-a) [[ 



C=b j£ V. (Ps - a) (p« + b)J fj-j \ {p s +a) (p s + b) 

for some functions c\ and c 2 independent of rife. Since the same argument may be carried out in all lattice 
directions it follows that ci and c 2 are constant. Using (|2.32[) . and the fact that lF(a,b) = p(a)zp(b)z, 
equation (|2.37a[) at £ = b then becomes 

(Pfc + &) ) U = ft« - (Pfc -&*)« + ci^(a, -6) (^j) + c 2 F{-a, -b) (^^) • (2-39) 

Writing 

U = AF(a, 6)S^ + SJ(a, -&)S£ + £F(-a, 6)§e + UF(-a, -&)§£, 
and taking the coefficient of JF(a, b) in (|2.39|) one obtains 



(Pk + b) 



A(&) \ = / (p fc + q)(p fc + &)Syi- (Pfc-& 2 )§A 
A(6) J I §^ 



A crucial observation is that one is free to interchange the roles of a and b in the N lax matrices (|2.3p . and 
thus we also have 

'A(o)\ = / (Pfc + ^fa+^^bg-a 2 )^ ^ 
A(a)i I • 



(Pk + a) 
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Combining these gives the two expressions 

,2 ^ f§A \ _ (Pk + a )A(a)A(fe) - (Pk + 6)A(o)A(6) 

[ } \MJ~ MajMb) (2 ' 40a) 

, 2 tfJ&A (Pk - &)A(q)A(b) - (p fc - o)A( tt )A(6) 

(a " M l^J = wm ■ (2 - 40b) 

One final Wronskian identity results from considering equation (|2.38|) for A(a) and A(6). By eliminating the 
middle term one obtains 

(p* - fc)A(o)A(6) - (p fc - o)A(o)A(6) = (p fc + o)X(o)A(6) - (p fc + 6)A(o)X(6), 
and combining this with ()2.40|) gives 

&A = ^— b ({Pk + a)A(a)A(b) - (p k + 6)A(o)A(6)) . 

By considering the points C — ~b, C = a an d C — — a one can obtain expressions for the remaining functions 
§s, Sg and §2) in a similar fashion. □ 

This concludes the derivation of the 1ST for Q3^. The entire procedure may be summarised as follows: 
Given a staircase T in an N-dimensional lattice of the form stipulated in Subsection 12 .21 one then specifies 
an initial profile i>i along this staircase. This potential must satisfy the summability condition (|2.22j) . By 
solving equation (|2.13|) one can construct the Jost solutions along V, and thereby determine the functions A 
and B and the constants Cj. Finally one can use the singular integral equation (|2.31|) to determine A as a 
function of all lattice variables n±, ...,nN and obtain the solution u of Q3^ using Theorem 12. 101 

3. Degenerations 

Theorem 12.101 allows one to express the solution (|2.34[) of Q3a in terms of A(±a) and A(±6). We now 
consider the problem of degenerating this solution to the remaining equations in the lists (jl.ljl and (|1.2I) . 
To do so we follow the degenerations given in }20) which are limits on the parameters a and b and the 
dependent variable u, in which a small parameter e is introduced, and all degenerations are obtain in the 
limit e — > 0. Further choices must be made for the degenerations of the constants appearing in the solution 
in order that the resulting solutions be of the required form. These choices are given in [20] and are 
repeated here for convenience. The most complicated is the Q3^ — > Q2 degeneration, in which we make the 
substitution b — a(l — 2e). The summability condition (|2.22[) on the potential v guarantees the existence 
and differentiability of A at £ = ±a and £ ± b, and after taking the required limit on b we have 

p r + a 
p r — a 

and thus the solution to Q2 will satisfy the summability condition 

\vi\X l i ] (l + \i\)<oc. (3.1) 

i— — oo 

This is also the case for the solutions of Qls and H3. The solutions to H2 and HI, which are obtained by 
degenerating the Q2 and Qls solution respectively, require the parameter a to approach infinity. By choosing 
DCi such that 3Ci(o) — > 1 as \a\ — > oo then the solutions obtained will satisfy the summability condition 

+oo 

J2 N(l + |*|)< oo. (3.2) 

i— ~ oo 

In the following degenerations we expand the Jost function A(£) using Taylor series. For example after 
setting b = a(l — 2e) in the Q3a — > Q2 degeneration one writes 

A(6) = A(o) - 2aA»e + 2a 2 A"(a)e 2 + ... . 

Note that in the following the notation A' (—a) denotes A'(£) = g|-A(C) evaluated at ( = —a. In addition to 
this one must make use of the Wronskian identity 

( Pk + o)X(o)A(-o) - (p k - o)A(a)A(-a) = 2a, (3.3) 



3C (a, b) — > %i(a) > max 



Pr 
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obtained by combining equation (|2.38[) for A(a) and A(— a), as well as the two further identities obtained by 
differentiating (|3.3[) once and twice with respect to a. Degenerations of the exponential terms J- are easily 
computed and can be seen in |20j . For the degenerations where the parameters a or & approach infinity we 
make use of (|2.25[) which gives 

A ~ 1 + e(A e ) + e 2 (A £ 2) + O (e 3 ) as e -> 0, 

where functions A e and A e 2 may be calculated by taking the limit |£| — > oo of A(£). For each degeneration we 
state the limits that one must take on the various terms and give the final degenerated solution. The various 
new functions obtained involve shifts in an arbitrary n^-dircction, denoted again by with associated 
parameter pk- 

3.1. Q2. 

As a first case we consider the degeneration to Q2 which lies just below Q3,5 in the ABS hierarchy. The 
degeneration is 

— (- 

4a 2 \e 

and the choices that we make for the constants in order that the solution degenerates in this manner are 



a(l-2e), u-> — ( - + 1 + (1 + 2u)e 



Note that the four constants A, 23, C and 2), constrained by (|2.8|) have been replaced by the three new 
constants A, D and £ D with no constraint. After this degeneration one finds the new inverse scattering 
solution of Q2 to be 

u=i((£ + £o) 2 + i) +a^ + ^o)s(-a,a) + a 2 Z{a)+AV (3.4) 
+ -Ap(a) (l - 2aS{a, a)) + -T>p(-a) (l + 2aS(-a, -a)) (3.5) 



where 



S(a, b) = 1 (l + ( Pk + a)A(o)A'(6) - (p fc + &)A(a)A'(&) - A(a)A(&)) 
Z(a) = ^([(Pk + o)A'(o) + A(o)]A'(-o) - [(p fc - o)A'(-o) + A(-a)] A'(a) 

iV 

£ = £(a) = 2a^ 



N 

Pr 



9 ? 



3.2. Ql 5 . 

The degeneration from Q2 to Ql« is given by 



S 2 1 
—u 



4e 2 e 

and the required degenerations of the constants A, D and £ are 

2A _ 2D ^ ^ 2S 
e e e 

One then finds the solution of Qlg to be 

u = 23 (*£ + Co + 2aS(-a, a)) + A/j(o) (l - 2aS(a, a)) + 2)p(-a) (l + 2aS(-a, -a)) , (3.6) 

where the constants A, 23, CD and £ are constrained by 

<5 2 

25 + AAD = — . 

4 
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3.3. H3,5. 

To obtain the H3a solution we degenerate from the Q3.5 solution (|2.34l) by setting 

1 3 fVs 



and 



b =72> u ^ e — u 



A ^ e sM A )^ v^ e 3M%) i e^ £ 3 [^eV 



2 J V 2 / V 2 / V 2 

The solution which emerges is 

u =(a+ a'B) tfA(et) + (aG + £>) tT 1 A(-a) (3.7) 

where 

and the constants A, 25, C and D are constrained by 

yiD - se = — . 

4a 



3.4. H2. 

In order to obtain the H2 solution we choose to degenerate from the Q2 solution (|3.4[) by setting 

1 1 2 

a = — , u — > — hew. 

e 4 

The degenerations of the constants A, D and £ are 

A^A(e + 2g 1 e 2 \ V -> .A (-e + 2 0l e 2 ) , £ -> 2eo , 
where .A, g and pi are unconstrained constants. The solution of H2 is then 

u=(g+g ) +2(g + g )A c + 2A £ 2 - A 2 + g2a(q + Oi+A e Y (3- 



where 

JV 
r=l 

3.5. HI. 

Finally we derive the inverse scattering solution to HI by degenerating from the QI5 solution (|3.6[) . Wc 

set 

a = — , it — > e^u, 
e 

and require that the constants degenerate to 

S-^S, X->|yt(l + 2eie), D->^(-H-2eie), £ G -> 2e£ OJ 

where 

yi 2 - S 2 = -1. 

The solution of HI which emerges is 

u = 'B(Q + Q a + A^ +aA{g + Q 1 +A e y (3.9) 

This solution agrees with that given in [6] however is formulated in terms of the expansion of A(£) for large 
£ rather than using the Gel'fand-Lcvitan approach. 
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4. Reflectionless Potentials 

Consider the case of an arbitrary reflectionless potential, i.e. one for which = 0, where A vanishes at 
exactly one point C = k in CP*. The integral equation for A(£) becomes 

where we have written C\ = 2kp° . The solution is 

i + a(Op°p(A0 



A(C) = 



1 + p°p{k) 

C+fe 

^ (a+k){b+k) 



where A(£) = frr, and thus by ([OS)) we have 



= S(o, 6) = 



1 i (a-k)(b-k) o (i\ 

1 + ta+kY(b+k)P py K ) 



1 + /9°p(/c) 

For the reflectionless case where A vanishes at M points k\, kM in J 1 *, one can show that the solution of 
(EOTD is 

A(C) = l-c T (l + M) Ui + k)~ t, 



where c T = \ Ci, C2, CWj , r = ^/)(fci), p(k]\i)j , ]fy = fcjJjj and the Cauchy matrix 3Vt is given by 

p(ki)Cj 



k% -\- kj 



This is in fact the quantity V(C) given in |20j . and from the various identities given in this paper one can 
show that 



§ A =§(a,b) = c J [al + Kj [I + Mj [bl + Kj r. 

This is exactly the N-soliton solution obtained by the Cauchy matrix approach in |20j . Since the degenerations 
follow those from this paper, any reflectionless potential will give rise to pure soliton solutions to all equations 
in the degeneration scheme, that is all ABS equations below Q4. 

5. Other Lattice Equations 

In this section we give a short discussion on how the machinery presented in the previous sections allows 
one to obtain inverse scattering solutions to some previously known lattice equations, including some which 
are not multidimensionally consistent. We begin by deriving a closed- form lattice equation for the quantity 
A(£). Consider the evolution equation (|2.5j) for A(£) taken along opposing sides of a lattice element in the 
n- and m-dircctions: 

( ff + C)fe + (P - 0^ = ^-^-f)u-Pu (5.1a) 
A(C) A(C) It 

, + fl |<O + (,- C )*K) = ^-^ (5.1b) 
A(C) A(C) U 

One can then show that if u solves Q3^ then the quantities on the right-hand side of these equation are in 
fact equafl, and thus equating these shows that A(£) satisfies 

,„ + c) |k) _ {P _ Q m = iq+Q m- ( ,- m. (u) 

'HO a (0 a(0 a(0 

Solutions to this equation exist whenever £ ^ CP, provided that the potential v satisfies (|2.22[) . Let a p :— 
(p 2 — C 2 )^; a q (p 2 ~ C 2 )'j then in terms of the Jost function ip this becomes 

" P (<p((W(o - v(cmo) = oi q (ipicwio - £(cmo) , (5.3) 



this is the (1,2) element of the consistency equation for the Lax pairs in the n- and m-directions 
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that is <f(C) satisfies the lattice potential modified KdV equation [23J [2S] with parameters a p and a q . By 
setting W(C) := ^jf} we then sec that W satisfies 

v(C) 

W{Q = (a p W(Q-a q )(a p -a q W(0) 
W(C) (a p W(() - a q ) (a p - a q W{()) ' 

which is the lattice modified KdV equation [22]. Equations (|5 .3[) and (|5.4p are invariant under the transfor- 
mation ( — > —( and thus we have inverse scattering solutions to these equations for all £, independent of the 
choices of the parameters a and b. By letting these parameters approach infinity the solutions obtained will 
give rise to potentials satisfying (|3.2|) . Of course equivalent equations to (|5.2[) . (|5.3[) and (|5.4|) hold in any 
pair of lattice directions, and each of these equations is multidimensionally consistent. If we now break the 
covariance between the n- and m- lattice directions in ()5.2|) by setting ( = p, we have the following lattice 
equation for the quantity A(p): 

2 P = ip + q M + i P - q) m. (5.5) 

Hp) A(p) 

The (weaker) multidimensional consistency properties of this equation arc discussed in [20] . By setting C, = p 

in (|5.1a[) and £ = q in (|5. lb|) . then equating the right-hand sides we have A(p)A(q) = A(q)A(p) which shows 
that it is self-consistent to express A(p) in terms of the r- function: 

A(p) = I A(q) = L 
Using this in (|5.5p gives the 6-point equation 

2pff=(p + q)£f + (p-q)fl 

By supplementing this with the similar equation obtained by setting f = q in (|5.2|) . one can obtain Hirota's 
discrete-time Toda equation [15] : 

(P - qfff- (p + q?IJ+ 4pqf = 0. 
We now derive two final lattice equation of KdV type. From the QI5 equation, if one makes the degeneration 

a = e, u —> eu, 



then the equation becomes 



i»-«)(;-g) = 4, (5.6) 

(u — u)(u — u) P 

which is the Schwarzian KdV equation [22], also known as the cross-ratio equation (or Q1 Q ). By choosing 
the constants such that 

23 H 23 

A^-, V^-, eo^l-2e M o 

with fj. a constant, the inverse scattering solution to (|5.6[) is then given by 



u = %(u + h -'T), (5.7) 



where 

M 



T=f(A(0)A"(0)-A(0)A"(0)) -A(0)A'(0), M = E ~ 



2 V / 1 Pr 

This solution exists whenever holds with 3Ci(0) 7^ 0. Alternatively from the solution (|3.9[) of HI, also 
known as the lattice potential KdV equation JTBJ [23] [2S], if one defines ui := u — u then the quantity w 
satisfies the lattice KdV equation [TB] 

UJ-^ = {p 2 -q 2 )(l-i\ (5.8) 
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to which we have the inverse scattering solution 

u = (S - (-l)"+"Vl) (p - q + A e - A £ ) , (5.9) 
which exists provided (|3.2[) holds. Note that equation (|5 .8[) is not multidimensionally consistent. 

6. Conclusion 

In this paper we have developed a discrete inverse scattering transform for a hierarchy of integrable lattice 
equations. At the top of this hierarchy lies Q3^, for which the discrete 1ST was carried out. The forward 
scattering problem was carried out along a multidimensional staircase within an N-dimensional lattice. Even 
though this staircase may traverse multiple dimensions, the forward scattering problem along the staircase is 
one-dimensional. Building on the machinery developed in [5J the multidimensional complex- valued solution 
(depending on N independent variables and N distinct lattice parameters) was shown to be expressible in 
terms of the quantity A(£), which is obtained by solving a singular integral equation which is a discrete 
analogue of the singular integral equation obtained in the 1ST for the continuous KdV equation [15] . It was 
also shown that the soliton solutions presented in |20| correspond precisely with the solutions arising from 
reflectionless potentials. Following the degeneration scheme given in [5D] the solution to Q3,5 was degenerated 
to give inverse scattering solutions to all remaining ABS equations (except for Q4, which will be addressed in 
due course), as well as some other previously known integrable lattice equations, including the lattice KdV 
equation which is not multidimensionally consistent. 

The summability condition (|2.22[) placed on the potential is the only requirement needed to carry out the 
discrete 1ST, provided that it is scaled so as to be admissible in accordance with Definition ^. 41 This condition 
changes throughout the degeneration scheme, and the condition placed on the Q3^ potential reduces to (|3.2[) 
for the cases of H2 and HI, which agrees with that obtained in the case of the continuous KdV partial 
differential equation, and the condition found in [B]. 
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